Simple equation of state for hard disks on the hyperbohc plane 
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A simple equation ol state for hard disks on the hyperbolic plane is proposed. It yields the exact 
second virial coefficient and contains a pole at the highest possible packing. A comparison with 
another very recent theoretical proposal and simulation data is presented. 



It is well known that hard-core systemsi represent use- 
ful models that allow both the derivation of some rigorous 
results in statistical mechanics as well as the computation 
of some particular quantities such as virial coefScients. 
For monocomponent hard-core systems there cannot be a 
gas-liquid transition due to the lack of an attractive part 
in the intermolecular potential, but they show crystalline 
and/or amorphous phases and the way these phases arise 
and even their actual existence are still open problems. 

Insight into the thermodynamic behavior of hard-core 
systems has in the past been sometimes gained by con- 
sidering similar systems in higher dimensionsi^ The ra- 
tionale is that one may obtain a rough idea of the ther- 
modynamic behavior of say three-dimensional hard-core 
fluids at high density by looking at models in higher di- 
mensions in which the same phenomenology is present 
but at a lower density and thus the problem may become 
mathematically more tractable. In a different context, 
but nevertheless dealing with a somewhat related ques- 
tion. Modes and Kamicn^i have addressed the connection 
between the equation of state and disordered packings of 
hard disks. Since in simulations of monodisperse hard 
disks it is difficult to prevent crystallization, the consid- 
eration of a curved surface, in this case a hyperbolic plane 
with a curvature near a known regular tesselation, serves 
to frustrate global crystalline order and allows them to 
find the equation of state via molecular dynamics. Fur- 
ther, they have also developed a free area theory for the 
packing derived from the nearby tesselation. Studies in 
the hyperbolic plane in connection with glass-forming liq- 
uids and bulk behavior of physical systems have also re- 
cently been reportedi^ 

A few years ago we proposed a simple and accurate 
equation of state for a hard-disk fluid.— This equation is 
built so as to yield the exact second virial coefficient and 
also to have a single pole singularity at the close-packing 
fraction, namely 

p/pknT =[l- b2T] - (1 - &2fyraax)fyV^max] ^ > (1) 

where p is the pressure, p is the number density, A:b is 
the Boltzmann constant, T is the absolute temperature, 
62 = 2 is the reduced second virial coefficient, 77 = ao{a)p 
is the packing fraction, with ao{a) = (7r/4)(T^ the area of 
a hard disk of diameter a, and r?max = V^tt/G ~ 0.9069 
is the value corresponding to crystalline close-packing. 



The major aim of this Note is to provide an answer to 
the question of whether a proposal such as Eq. ([1]) , prop- 
erly generalized, may also be useful for hard disks on the 
hyperbolic plane. 

The development goes as follows. First, we note that 
on a two-dimensional manifold with constant intrinsic 
curvature K < 0, the area of a disk of diameter a is 
aK(cr) 27r|A'|-i [cosh(|A'|i/2cr/2) - 1] and that the 
packing fraction is given by rj — aKio')p. Further- 
more, on such a manifold the associated (reduced) sec- 
ond virial coefficient depends on the (reduced) diameter 
|A:|i/2cr, namely b2{\K\^/^(j) = aK(2cr)/2aif(cr)4 Thus, 
for a given value of the reduced diameter \K\^/^a, the 
only other requirement in our formulation is the corre- 
sponding value of rymax, irrespective of the fact that the 
resulting configuration be an ordered one or not. As 
far as we know, this value is only known analytically 
for the so-called {p, q} tesselations (restricted in the hy- 
perbolic plane by the condition 1/p + 1/q < 1/2), in 
which the packing corresponding to the highest possible 
density and the associated reduced diameter are given 
by ?7max = [cos(7r/g)/sin(7r/p) - 1] /(p/2 - 1 - p/q) and 
|A:|i/2cr ^ 2cosh"^ [cos ( tt/ g) /sin (tt/p) ] , respectively.— 

The top panel of Fig.[T]shows the compressibility factor 
Z = p/ pksT as a function of the packing fraction t], ob- 
tained from the free area theory of Modes and Kamien^ 
and from the use of Eq. ([T]) , in the case of three isostatic 
tesselations, that is, tesselations with p = 4. Clearly, 
the compressibility factors in both approaches are quite 
close for high packing fractions, but they differ substan- 
tially at low ones. We also compare the results for the 
packing fraction dependence of Z obtained with (to our 
knowledge) the only presently available simulation data 
of Modes and KamienS for two values {\K\^^^a = 1.062 
and \K\^^^a = 1.060) of the hard-disk diameter which 
are close to the one corresponding to the isostatic {4, 5} 
tesselation, namely \K\^/^a — 1.0613. The performance 
of both theories with respect to the simulation is rather 
poor with, if any, a very slight superiority of our ap- 
proach. Nevertheless, no definite statement about such a 
comparison should be made at this stage given the very 
small number of particles (1 to 9) used in the simulations 
and the intrinsic difficulties associated with simulations 
on the hyperbolic plane.— 

To further test the performance of both theories in a 
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FIG. 1: (Color online) Top panel: Compressibility factor 
Z as a. function of 77 for the isostatic {4, 10*^}, {4, 10}, and 
{4,5} tesselations {\K\^'^<J = 1.7628, 1.6169, and 1.0613, 
and ?7max = 0.4142, 0.5750, and 0.7206, respectively); solid 
lines: Eq. lU; dashed lines: free area theory of Ref. H the 
symbols are simulation results also obtained in Ref. y for 
\K\^/'^a = 1.060 (squares) and \K\^/^a = 1.062 (diamonds). 
Bottom panel: Compressibility factor Z as a function of 77 
for the (Euclidean) hexagonal {6, 3} tesselation (|if|^''^cr = 0, 
'7max = \/37r/6); solid line: Eq. (Il|; dashed line: free area 
theory of Ref. |3|; the symbols are the simulation results from 
Ref. 0. 

more controllable situation, in the bottom panel of Fig. 
[T]we present the results obtained in the Euclidean limit 
{\K\^/'^a = 0) using the hexagonal {6,3} tesselation. Of 
course in this limit in our approach we recover the results 
of our earlier proposal for the equation of state of a hard- 
disk fluid ^ which we had proved to be rather accurate. 
The free area theory;^ on the other hand, performs very 
poorly. 

While the philosophy behind the derivation of our 
equation of state is totally different from that of the free 
area theory of Modes and Kamien^^ they both share the 
property of having a pole at 77 = ?/max- In fact, from 
Eq. ^ it follows that Z ~ [2 - 62(|A'|^/2cr)?y,„ax]~^(l - 
??/??max)~"^ while the free area theory yields Zfteo area — 
2(1 — ry/rymax)"^ for 77 close to Tymax- Note that if 



lifji/V = 1.0613 then [2 - b2i\K\^/'^a)r]„^]-^ ~ 2.20. 
This explains why at high packing fractions the numer- 
ical values of the compressibility factor obtained with 
both approaches are not all that different. One should 
point out that the free area theory has been constructed 
with the particular aim of performing well at high densi- 
ties while ours attempts to capture both the low and high 
density limits. On the other hand, Zf^^c area is nonana- 
lytic at 7^ = (for the above value of \K\-^/^a, Zfrco area = 

I.26277I/2 4. 1.432877 1.594177^/2 H ) and incorrectly 

predicts Zfy-cc arca(O) = 0, whereas by construction Z 
given by Eq. ([1]) yields the exact second virial coefficient 
b2i\K\^^^a) and also allows us to estimate the higher or- 
der virial coefficients bni\K\^^^a). Thus, from Eq. ([T|) 
and again for \K\^^^a = 1.0613, we get 63(1.0613) = 3.55, 
64(1.0613) = 5.36, and 65(1.0613) = 7.76. Compar- 
ing with the numerical values obtained by Modes and 
Kamien^ for \K\^/^a = 1.10, namely 63(1.10) = 3.39, 
64(1.10) = 4.62, and 65(1.10) = 5.83, we find that our es- 
timates are of the right order of magnitude although their 
accuracy seems to be worse than the one we got for the 
Euclidean casci^ As a final comment one should add that 
Eq. IH) may lead to (unphysical) negative values of Z in 
the interval 77max/ [b2{\K\^^'^a)ri^a,^ - l] < 77 < 77inax if 
b2i\K\^^^a)T]niax > 2. This may well be a limitation of 
our formulation but the point certainly deserves further 
study. 

In summary, in this Note we have presented an exten- 
sion of our former equation of state for a hard-disk fiuid 
to deal with the same system on the hyperbolic plane. 
In contrast with the free area theory of Ref. d, which is 
rather more complicated and devised to perform well in 
the high density region only, the main assets of our pro- 
posal are its simplicity and the fact that it caters both for 
low and high packing fractions. Our results indicate that 
it might be as accurate and even improve on the perfor- 
mance of the free area theory of Modes and Kamien.— A 
deeper assessment of its full value is not possible at this 
stage due to the lack of simulation results for different 
values of \K\^^-^a. We hope that our simple theory will 
encourage the performance of such simulations. 
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^ For a very recent review of the properties of these systems 
one may look at Theory and Simulation of Hard-Sphere Flu- 
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ids and Related Systems, edited by A. Mulero (Springer, 
Berlin, 2008). 
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